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^3 ■ Abstract 

We investigate what can be learnt about the rj -rj' mixing by means of dispersive rep- 
resentation of axial anomaly. We show that our method leads to the strong bounds for the 
rj — rf mixing angle: 9 = —15.3° ± 1°. Moreover, our result manifests also a dramatic 
dependence of the width T^-^ on the mixing angle 9. This property explains how the 
relatively small mixing strongly effects the decay width. 



1 Introduction 



oo 

The problem of the t]T]' mixing and particularly the value of the mixing angle 9 has attracted a 
lot of interest for a long period of time. In the last decade the experimental data significantly 
improved (see e.g. 0j and references therein). So the improvement of the accuracy of the theo- 
retical prediction of the mixing angle 9 seems to be important. Very rich literature is dedicated 
to this question (see ED- 031). 

The well-known estimation, based on Gell-Mann-Okubo mass formulas give the value of 
9 about -10° for a quadratic mass formula and about -23° for a linear case. Obviously, the 
accuracy of these estimations is not high. 

In the well-known paper [2J many years ago by use of the expansion of the composite 
operators in the interpolating fields the values of the 777/ — > 27 and J/ij) — > i]^')^ were con- 
sidered in the terms of the octet currents. The value of the mixing angle 9 was obtained there 
by comparison of theoretical results with experimental results (existed at that time). Two pos- 
sible mechanisms for the pseudoscalar field mixing, offered in J3), so called mass mixing and 
current-mixing models were discussed, and the values of 9 about — 15°± 1.8° and — 19.2°±2.2 
correspondingly were obtained 
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The analysis of the axial anomaly generated decays r](r)') — > 27 was also performed in fl4) 
(in the framework of ChPT) and [|5]| and the estimate 9 = —20° -j- 25° was obtained. The compe- 
hensive general analysis of the pseudoscalar mesons decays by making use of both abelian and 
nonabelian anomalies was performed in (8]|. This approach is based on a standard (differential) 
form of the anomaly and on the pole approximation. 

Another systematic anomaly based approach was developed in (6JI3, where a large set of 
decay processes was investigated. The mixing angle was estimated as 9 = — 17° ± 2°, which 
significantly differs from the estimation of [H: 9 = — 9° . The possible explanation of such 
discrepancy is a different role played by the continuum contributions in this two approaches. 
The significance of the continuum contribution was explicitly mentioned in 0. It was noted, 
that for T]T]' — > 27 decay the value 9 = — 17° ± 2° can be obtained only if one takes into account 
the continuum contribution which cancel the effects of the large SU(3) breaking (i.e. the ratio 
f&l ' f-w = 1-25), otherwise the value of 9 grows up to -21°. 

The effect of isoscalar pseudoscalar continuum on the mixing angle was also examined in 
[fTTTl . In this approach the masses of the isoscalar mesons and their decays were discussed. 
The value of the mixing angle was estimated in the range —30.5° < 9 < —18.5°. The value 
—17° < 9 < —13° was obtained in ifTOl . where phenomenological analysis of a set of decays 
was performed. 

All the above results were obtained in the mixing scheme with one angle. The case of the 
two mixing angles tfT6l [T2l [T3l [T4| we will discuss in the section 4. 

We can see that all results are compatible , but with relatively large uncertainties. 

That's why it seems reasonable, and this is the main aim of this paper, to use the dispersion 
form of axial anomaly to find some model-independent and precise restriction on the value of 
the mixing angle. We will develop the method, used in [|23l . where a very precise prediction on 
7T° — > 27 decay was obtained. 

In section 3 we will discuss the one angle octet-singlet mixing scheme and a very precise 
restrictions will be obtained. The scheme with two angles will be discussed in sect. 4. 



2 77 —> 27 Decay Width Discrepancy 

Consider a matrix element for a transition of the 8th component of the axial current into two 
photons with momenta p, p' and polarizations e a , e^: 

T pa p(p,p') = (p,e a ;p',e' p \jj%\0) , (1) 

where 

j{ $ = -^(«7m75W + hiflhd ~ 2s7 M 7ss) • (2) 

Here u, d and s are fields of u,d and s respective quarks. 
We are interested in the decay into two real photons, thus the general form of T Ma/ g(p, p') can be 
expressed as: 

Tp a p(p, p') = Fx (q 2 )qp€ a p pa p p p' a + ^F 2 (q 2 )(p a e^ pa - p' p ep apa )p p p' a , (3) 
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where q = p + p'. 

Using the dispersive approach ifTTll (see also [12711 and references therein) one can derive the 
exact "anomaly sum rule" lfT8l [T9l |20ll : 



(4) 



oo 

J ImF 1 (q 2 )dq 2 = V2a(e 2 u + e 2 d -2e 2 s )N c =^c 



o 



where e u = 2/3, e d = e s = 1/3 are electric charges of u, d, and s quarks, N c = 3 (the number 
of colors). 

Notice, that in QCD this equation does't have any perturbative corrections 11221 . and it is ex- 
pected that it does not have any nonperturbative corrections as well due to t'Hooft's consistency 
principle [Tl9l . It is also important that at q 2 — > oo the function ImFi(q 2 ) decreases as l/q A . 

Based on the dispersion relation © for an axial anomaly, a high precision prediction for 
7T° — ► 27 decay was obtained ll23l ( 1.5% , mixing 7r° - 77 was taken into account) and later was 
confirmed by experiment. 

Let us first try to saturate the above relation by the 77 contribution only. We use the following 
standard definition of the 77 decay constant: 

(0\jS\v)=if^- (5) 
The general form of the ^-contribution to T^^p, p') can be written as: 

T^(p,p') = ~f v ~2Z — ^q^afiXoVxV'c ■ (6) 
q 

where A v is a constant. So we see, that in the naive approximation, when only 77 contribution is 
accounted in the l.h.s. of the sum rule for anomaly relation (HJ), one can find A n 

A 1 

Then one can easily calculate the decay width rj — > 27 ETI : 

~ 1 a 2 

v ^< = 3 32^ ~f* ■ (7) 

If we substitute experimental numbers of a, m v and f v = 1.2 f n ~ 150MeV we'll get the 
numerical value 

r^2 7 = 0.12 keV , 
which is in a striking contradiction with an experimental value 

r^ 27 = 0.510 ±0.026 keV . 

Such a large contradiction with exact anomaly dispersion relation motivates us to consider 
the effects arising from the other states contributions to the sum rule. We will mainly discuss 
the mixing of 77 and 7/ mesons, because calculation show, that effect of n — 77 mixing is found 
to be extremely small. We will say a few words about this later. 
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3 One Angle 77 — 77' Mixing Scheme 

Let us now consider the decay 77 — ► 27 using the dispersive approach and taking into account 
mixing of 77 and 7/ mesons. We start with the well-known octet-singlet mixing scheme. 

Following E4ll25l . let's introduce nonorthogonal states |P 8 ) and \P ) and the corresponding 
fields ips, tpo, coupled to jfg and 

(0\J$\Pi) = i5 kl f k qr, fc=8,0. 

Nonorthogonality of the fields y> 8 , <^ corresponds to the non-diagonal term AH = mL,</? 8 </? in 
the effective interaction Hamiltonian. In the presence of such term the standard PCAC relation 
is modified in the following way: 

d n J £ = fv( m2 v ¥8 + m l v 'Vo) , (8) 
The fields y2 8 ,yj are expressed through the physical fields ip v , </y as 



ip 8 = tp v cos 9 + tp v > sin 9 , (9) 
<fo = —(p v sin 9 + <f V ' cos 9 . (10) 



Mixing angle 9 can be expressed from (|9l).(fT0l) in terms of masses as: 

2m 2 , 

tan20 = „ m - . (11) 

Now ImF 1 (q 2 ) is given by the sum of contributions of 77 and 7/ mesons. In order to separate 
the formfactor Fi(q 2 ), multiply T fia p(p, p') by q^/q 2 . Then, taking the imagenary part of Fi(q 2 ), 
we get: 

ImF 1 (q 2 )=Imq fl ^(2 1 \J ( l S | 0) 

= — |/m(27 I m 2 (Lp v cos9 + (p v i sm9) + m 2 vl (—cp v sia9 + cp v i cos6>) | 0) 

2 

Tn 

= Trf v [A v cos95(q 2 - m 2 ) + A v > — |- sin6'5(g 2 - m 2 ,)) 

m r 1 ' 

2 2 

- A^-sm95{q 2 - m 2 ) + ^^cos^(g 2 - m 2 )] , (12) 



where A v ( A n t) is the amplitude of decay 77(77') — > 27. 
If we employ sum rule ©, we obtain: 



2 2 2 

m n . „ . m rm' . „ . m W „, /2 



tt/„[A,cos0 + A,'— f sin0- A,— ^f-sinfl + A^— f-cos0] = W -a . (13) 



m~, m~, 



4 



Now let's express amplitudes in terms of decay widths, employ the relation (TTTT) for a mixing 
parameter nil , and finally get the equation for a mixing angle 9: 



ml 1 , ml, 8 ml 



cos# + /?^sin#- -(—I- - 1) tan 20 sin # + -(1 f ) tan 20 cos# = £ , (14) 

m^, 2 m^ 2 m^, 

where the dimensionless parameters (3 and £ were introduced: 




^ = M {= / ^ i _ (15) 



V 967T 3 r^ 27 /2 



Before starting a numerical analysis, let us discuss the accuracy of the obtained equation 
and possible additional contributions. The other (besides 77 and rf) pseudoscalar meson states 
contributions can arise as some additional terms in the r.h.s of eq.(fl~2l). But due to the fact, 
that ImF\ decrease as l/q 4 at large q 2 , from dimensional ground one can easily conclude, that 
higher resonances contribution should be suppressed significantly stronger than (m v / 'm res ) 2 , 
more probably as a fourth power of the mass ratio. This fact effectively suppresses uncontrolled 
contributions of higher resonances and continuum and, consequently, allow us to reach a higher 
accuracy. This is one of the advantages of this method. The numerical analysis show, that 
possible influence of this uncontrolled effects on the solution of the eq. (fT4l) (i.e. mixing angle 
9) is very small (at a few per cent level). 

We also calculated the contribution of the n — rj mixing effects in the same way (in fact all 
what is necessary for this was done in [|23l ). The pion mixing contribution was found to be 
extremely small (less than a per cent, suppressed as (m w /m ri ) 4: ), that's why we don't mention 
it. 

For a numerical analysis of our anomaly sum rule equation (fl4l) we use the following exper- 
imental data (PDG Review 2008) [26J: m v = 547.853 ± 0.024 MeV, rry = 957.78±0.24MeV, 
r^27 = 0.510±0.026keV, T^^^ = 4.30±0.15 keV. Also for a convenience we will present 
in units of a decay constant of pion f n = 130.7 MeV. The only parameter in eq. (fT4l) that 
has no a well-defined experimental value, is f^. If we use a well established prediction of ChPT 
(see e.g. [QJ) /„ = 1.28^, we get: 



= -15.3°±l°. (16) 

Error bars include both theoretical and experimental uncertainties which happened of the 
same order. It is useful to express f v from eq. (fl4l) as a function of 9. Remarkably, the function 
f v (9) within the interval of reasonable values of shows quite strong dependence on the mixing 
angle 9 (fig.l). Moreover, expression (fl~4l) manifests also a dramatic dependence of the width 
F v — > 27 on the mixing angle 9. This property of the anomaly motivated equations (fT2l) . (fl~4l) 
explains how the relatively small mixing strongly effects the decay width. 

It is interesting also to note, that eq. (fl"4l) as an equation for 9 formally has 6 roots for 
fv ~ 0-5/tt an d only 4 roots for f v < 0.5/^ (see fig. 2). Moreover, the only physical root is 
among these two disappearing roots. 
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4 Two Angle 77 — 77' Mixing Scheme 



An approach with one mixing angle in study of rj — rf system dominated for decades. But in the 
recent years there was a rise of interest to the mixing scheme with two angles. It was noticed 
that taking into account the chiral anomaly through perturbative expansion in ChPT can lead to 
the introduction of two mixing angles in description of the 77- 77' system |fl6l|. There were also 
performed some phenomenological analysis of various decay processes in the approach of two 
angle mixing scheme |[T2l[T4l . 

In this section we will consider a two angle mixing scheme for the 77 — 77' system using 
the approach developed in the previous section. We suggest a straightforward generalization of 
approach ll24l l25l by combining PCAC relation ([8]) with presence of two mixing angles 9\ and 
Mcf. El): 



(p 8 — (p v cos 9 2 + Vti' sin 9\ , (17) 
cpo = — ip v sm9 2 + <Prf cos#i . (18) 

When mixing angles 6>i and 9 2 are equal to each other, one comes back to the scheme with one 
mixing angle. The anomaly sum rule equation for a two angle mixing scheme now reads: 

22 2 

cos 9 2 + 0^1 sin 9 X - sin °2 + P^f cos 9 t = £ , (19) 

m*, m~, 



where (3 and £ are defined by (fT5l) and 



m 2 , sin 26*1 — nil sin 29 2 

^2 n x v 

2 (cos 9i cos 9 2 — sin 9\ sin 9 2 ) 



m z , = — -2 . (20) 



If we use the experimental values for masses m v , m r f and widths T,^^, T v /^ 2l as in the 
previous section, we will get the equation which contains as parameters angles 9\, 9 2 and f v . 
For physically interesting values of angles we get the plot shown on the fig. 3. The resulting plot 
for a full space of parameters is presented on the fig. 4. We may see that the angles play rather 
different role. While 9\ is tightly bounded by the data, like in the case of one angle mixing 
scheme, the sensitivity ro the changes of 9 2 is rather weak. 



5 Conclusions 

It was shown that the use of dispersive representation of electromagnetic axial anomaly leads to 
tight bounds for the 77 — 77' mixing angle: 9 = —15.3° ± 1°. This is model independent precise 
prediction which should be taken into account in future analyses. 

One of the main advantages of this method is effective suppression of the higher contribution 
and higher resonances which are usually the main source of uncertainties. 

Moreover, our result manifests also a dramatic dependence of the width T v ^ 27 on the mixing 
angle 9. This property explains how the relatively small mixing strongly effects the decay width. 

We considered also two angle mixing scheme and and found that rather weakly correlated: 
while the first angle is still tightly bounded, the second one may vary in a rather wide region. 
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Figure 1 : Mixing angle 6 as a function of decay constant f v in the one angle mixing scheme. 
Dashed lines correspond to the errors of the experimental data input. Dot-dashed horizontal line 
indicates the f v = 1.28/V level. 



<*? K i.o 




Figure 2: Mixing angle 6 as a function of decay constant f v in the one angle mixing scheme - 
the full parameter space. Dashed vertical lines indicate asymptotes. Dot-dashed horizontal line 
indicates the /„ = 1.28 f n level. 
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Figure 4: Two angle mixing scheme - the full parameter space. Different contours correspond 
to different f v (in units of f^, denoted by numbers). The dotted bisector corresponds to the 
reduction of the one angle mixing scheme (6i = 6 2 ) 
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